Bound Molecules in an Interacting Quantum Walk 
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We investigate a system of two atoms in an optical lattice, performing a quantum walk by state- 
dependent shift operations and a coin operation acting on the internal states. The atoms interact, e.g., 
by cold collisions, whenever they are in the same potential well of the lattice. Under such conditions they 
typically develop a bound state, so that the two atoms effectively perform a quantum walk together, rarely 
moving further from each other than a few lattice sites. The theoretical analysis is based on a theory of 
quantum walks with a point defect, applied to the difference variable. We also discuss the feasibility of an 
experimental realization in existing quantum walk experiments. 



PACS numbers: 42.50.Ex, 31.15.aq, 03.65.Ge 

Quantum walks are the discrete-time and discrete- space 
analog of free quantum particle motion. They exhibit many 
of the well-known features of quantum dynamics: In the 
translation invariant case their motion is governed by a 
dispersion relation, from which the position distribution 
at large times is obtained as the distribution of group ve- 
locity. By modifying the coin operation of the quantum 
walk on a small number of sites, one can create the ana- 
log of potentials, with non-trivial scattering cross sections 
and bound states. In this paper we consider the possibil- 
ity of an interaction between two walking particles, which 
becomes effective whenever the particles occupy the same 
site. In the simplest case, the Hadamard quantum walk, the 
interaction consists of just an additional phase. We show 
that with such interaction one can expect the formation of 
molecules, with all the features known from ordinary quan- 
tum dynamics of interacting two-particle systems: The two 
atoms stay together forever, performing a quantum walk in 
its own right. The distance between the particles remains 
bounded in the sense that they occupy nearby lattice sites 
with high probability, with estimates valid for all times. 

Quantum walks have recently been demonstrated exper- 
imentally in setups based on controlled neutral atoms [lj, 
ions EH, individual photons [4], wave guides El [6) or 
light pulses [7j. While all these experiments show the ex- 
pected one-particle dynamics, their ability to demonstrate 
the type of interaction we describe differs very much. In 
some realizations interaction is very natural if not unavoid- 
able: when two atoms sit in the same potential well of 
an optical lattice they are expected to pick up a collision 
phase depending on the duration of the contact. In spite of 
the fact that the atom-atom interaction would be much too 
weak for chemical binding in free space, the lattice creates 
a complex interference effect, which nevertheless leads to 
binding. The relevant conditions and parameters for an im- 
plementation of molecules of two atoms walking in a ID 
optical lattice are discussed below. In a suitable system 
molecules are easy to prepare by starting two particles on 
the same site. The initial wave function then typically has 
a good overlap with the bound state, and the group velocity 



of unbound pairs is higher than that of molecules. Hence 
one can just wait for the unbound atom pairs to move away. 
Our conclusion is that an experiment should be possible in 
the near future, perhaps even by extending an existing ex- 
periment. 

Our present work was partly inspired by analogous work 
in the continuous-time case IH, in which the motion be- 
tween neighboring sites is given by tunneling. The spec- 
tral structure of continuous-time lattice dynamics is indeed 
analogous to that of discrete-time walks, except for the fact 
that in discrete time the spectrum of the unitary one-step 
operator is on the unit circle, and in the continuous-time 
case the spectrum of the Hamiltonian is on the real line, so 
it makes sense to speak of high energies and low energies. 
In contrast, for discrete-time quantum dynamics there is 
no distinction between repulsive and attractive interaction 
potentials. Hence, in contrast to El, we cannot speak of 
"repulsively bound" molecules. In the discrete time case, a 
recent study of correlations in two-particle quantum walks 
|9j , numerically found enhanced correlations related to the 
molecule states studied in this paper. 

Free Quantum Walks. — Let us begin by describing quan- 
tum walks of a single particle. A walk runs on an s- 
dimensional lattice, which we take as the cubic lattice 
If . The walking particles generally have an additional 
quantum degree of freedom, described in a d-dimensional 
Hilbert space C d . We can think of this as the internal states, 
or a "coin" degree of freedom, which can be "flipped" by 
a suitable unitary operation. The basis states of the system 
Hilbert space are thus \x, a) with position x E If and in- 
ternal state label a. The dynamical time step is given by 
a unitary operator W, which commutes with lattice trans- 
lations. A crucial condition is locality: W has non-zero 
matrix elements only between states on lattice sites closer 
together than some fixed finite distance, called the neigh- 
borhood size of the walk. Since W commutes with trans- 
lation it is convenient to partially diagonalize this operator 
by Fourier transform, i.e., to look at the walk in momen- 
tum space. In the momentum representation wave func- 
tions depend on p E [— 7T, 7t] s , i.e., the Brillouin zone of 
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the lattice, with ip(p) E C d . The walk in this represen- 
tation acts by multiplying the wave function at p with a 
unitary d x d-matrix W(p). The locality condition then 
becomes the statement that each entry must be a polyno- 
mial in the variables e ±Wk , k = 1, . . . , s. The standard 
example is the so-called Hadamard walk in one dimension 
(5 = 1), for which 

w H { P )=s{ P) c=[ e ; e %)^{\ \) a) 

Here the first factor is usually called the (internal state- 
dependent) shift and the second the coin. In one lattice 
dimension a decomposition into (usually several) such fac- 
tors is always possible, but in higher dimension this is an 
open question. There is no intrinsic connection between 
the space dimension s and the internal state space dimen- 
sion d. In particular, it is not necessary to choose a different 
coin for each dimension (which would give d — 2 s [10]). 

For a non-interacting walk of two distinguishable parti- 
cles, we just take W 2 (pi,p 2 ) = Wi(pi) ® Wi(p 2 ) where 
Pi and p 2 are the momentum coordinates of the two parti- 
cles. Apart from a doubled lattice dimension (s f = 2s) and 
a squared dimension of the internal state space (d! — rf 2 ), 
this is a walk exactly as described above. Obviously, since 
W\ ® W\ commutes with the particle permutation, we may 
restrict the two particle walk to either the Bose or Fermi 
subspace. The influence of entangled initial states on the 
position distribution of the non-interacting walk has been 
studied theoretically in ifTTll . and experimentally with pho- 
tons in 0. 

Interacting Quantum Walks . — To make a two-particle 
walk interacting we need to introduce a dependence on 
the particle coordinates. The general way to introduce 
such modifications is to define a new transition operator as 
W — W 2 C, where W 2 is the unperturbed, fully transla- 
tion invariant two particle walk, and C is a space dependent 
coin operation, i.e., a unitary which acts (in a possibly dif- 
ferent way) on the internal degree of freedom at each site. 
We choose an operator C, which differs from the identity 
operator 1 only when the two particles are on the same lat- 
tice site (xi = x 2 ), in which case it acts on the joint coin 
space as some fixed operator T. We denote by N the pro- 
jection operator onto the set of collision points, i.e., the di- 
agonal in the x i5 x 2 -diagram. The overall interacting walk 
operator is 

W r = (W 1 ® Wi)((I — N) + FN). (2) 

Note that N only acts on the translation degrees of freedom 
and T acts on the coin space, so these operators commute 
and the second factor is unitary. Since T is taken to be the 
same for every collision point, we thus break the separate 
translation invariance, but retain the invariance of the walk 
under joint translations. Therefore, the total momentum 
Pi + P2 is conserved by Wr- 

We will develop the basic theory for such walks quite 
generally, for any T in any coin dimension d and also any 



lattice dimension (see Appendix A for the salient condi- 
tions). However, for illustration we will focus especially 
on the simplest case, which we call a walk with singlet 
collisions. This is defined by having the coin space at the 
collision points one-dimensional, so the coin acts just by 
a phase T = 7 1. In the Fermi case, and when the single 
particle coin space is two dimensional as for the Hadamard 
walk, this condition is automatically fulfilled even for the 
free walk. In the Bose case of the Hadamard walk, particles 
can only end up on the same site if they come from opposite 
directions, i.e., in a "symmetric singlet" ( | fi) + 1 it) ) / V2. 
If the combination of free coin and collision coin leave this 
state invariant, we are again in the singlet collision case. 

In order to define the interaction for more than two par- 
ticles, or for a quantum lattice gas of such particles, more 
work is required. We describe this in the supplement Ap- 
pendix B, mostly for the case of singlet collisions. 

Numerical Phenomenology. — We start the walk ([2]) from 
the initial state in which both particles are at the origin, and 
their joint internal state is antisymmetric. Since the Fermi 
subspace is invariant under W r , all subsequent collisions 
will also be of singlet type. Then by multiplying the initial 
vector t times with ([2]), and taking the modulus square with 
respect to the internal states, we get the joint position dis- 
tribution of the two particles at time t. The result is shown 
in Fig. [T] once for the non-interacting case (7 = 1), and 
once for the interacting case with value 7 = —1. Note 
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FIG. 1 . Joint probability distribution of two particles after t = 50 
steps of a Hadamard walk (antisymmetric subspace). Left panel: 
without interaction (7 = 1); Right panel: interaction with colli- 
sion phase 7 = — 1. The embedded (green) squares are for com- 
parison with the theoretically calculated peak velocity of the free 
Hadamard walk (left panel) and the peak velocity of the effective 
walk of molecules after Eq. ([9}. 

that due to Fermi symmetry both diagrams are symmet- 
ric with respect to reflection at the diagonal. The peaks 
near x\ — —x 2 are expected from the theory of the free 
Hadamard walk |[T2l . and are near ±y/2. The striking con- 
centration of probability near the diagonal is the hallmark 
of the bound state. The width of this distribution in the off- 
diagonal direction remains constant in time, whereas along 
the diagonal it shows the characteristic behavior of a ID 
walk. This will be proved below, and the peak velocity of 
this walk will be determined via Eq. ([9]) to be 1/3. For 
comparison, squares with edges at ±l/\/2 and ±1/3 are 
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drawn in Fig. [I] 

Turning from space and time diagrams to momentum 
and energy, let us consider the spectral properties of the 
operator ([2]). Since it commutes with joint translations, we 
can diagonalize it together with total momentum p. For 
each value of p we get an operator Wr(p), whose eigen- 
values we can compute. To treat this numerically, we close 
the system to a ring, which also discretizes p. The eigen- 
values of Wr{p), as a function of p, are shown in Fig. [2} 
Again, one part of this diagram is expected from the non- 
interacting case. Indeed, in that case W r (p) = W 2 (p) 
also commutes with translations in the difference variable 
(xi — x 2 ), which implies absolutely continuous spectrum 
for Wr(p)> These bands are seen in discretized form in 
Fig. [2j and the indication of continuous bands is that the 
spectral density increases with the size of the ring. In addi- 
tion, however, we see a single line in the gap between the 
bands. This is the bound state. 
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FIG. 2. Spectrum of the interacting walk operator for two par- 
ticles on a ring (28 sites). Left Panel: The non-interacting case: 
The relative momentum is then a good quantum number, and is 
represented as a third direction, orthogonal to the paper plane. 
For large rings the points fill the indicated surfaces. Right Panel: 
The interacting case (7 = — 1). The bands are very nearly equal to 
the projection of the left-panel figure onto the (p, uS) -plane. The 
additional feature introduced by interaction is the chain of simple 
eigenvalues in the band gap, i.e., the molecule states. 

Analysis: quantum walks w ith a point defect. — We first 
rewrite the problem in center of mass and relative coor- 
dinates. In momentum space this corresponds to the pas- 
sage from the variables (pi,p 2 ) to p — pi + P2 and 
k = (p 1 — p 2 )/2. Since p is conserved we can fix p and 
study the operator W r (p) as an operator acting on the k co- 
ordinate and the internal states. In this problem (for which 
p is just an external parameter) the interaction appears as a 
modification of a free walk acting only at the origin. The 
analysis of point defects in a general quantum walk is of 
some interest in its own right. We therefore drop the pa- 
rameter p for the moment, and consider this general prob- 
lem. So let 

W(k) = W 2 (p, k) = W x (| + k) ® W x (| - k) (3) 

denote the free walk matrix in relative momentum space, 
which corresponds, for each value of p, to a single particle 



ID walk with a four-dimensional coin space. We thus have 
to analyze the family of ID quantum walks Wr = W 2 (1 — 
N + NT) , where N is the projection onto the subspace at 
the origin, and T is the coin modification at the origin. 

Since N is a finite rank projection, the bands, i.e., the 
continuous spectrum of Wr, are the same as those of 
W 2 UHl Thm.IV.5.35]. To find the eigenvalues, suppose 
that ^(k) is the momentum space representation of an 
eigenvector with eigenvalue z. Then W(k)(l + (T — 
l)N)V(k) = z^(k). Note that is just the value of * 
at the origin in position representation, so N^(k) — i\) is 
a function independent of k. This ip determines the whole 
function \I/ by 

tf(fc) = (W(k) - z^WWil - T)^. (4) 

The condition for an eigenvalue z is thus that \I/ is normal- 
izable and ifj = 7V\I/. Since N is the projection onto the 
zero Fourier component, i.e., integration with respect to k, 
we introduce the operator 

and get the eigenvalue condition for ip in ([4]) to determine 
the eigenvector for eigenvalue z in the form R(z)(l — 
= This can be rewritten as 

r^= (l-R(z)- 1 )^. (6) 

Now it turns out that for z not in the spectrum of W, the 
operator on the right is unitary (see Appendix C). There- 
fore, the eigenvalue condition ([6]) can be satisfied for any z 
in the gap in Fig.[2j for a suitable collision operator T. 

Applying this to our standard example, the Hadamard 
walk with singlet collisions, the space N is one- 
dimensional, and the integral ([5]) is readily evaluated by the 
residue theorem, turning ([6]) into a formula for the collision 
parameter g in T = 7I = e ig l, given the total momentum 
p and the phase z = e lUJ . It turns out that this relation can 
be solved for u, giving 



e luJ = — _ - {^cosp±i^smp + 4(l - cos g)j 

provided sincj • sin(g — uj) > . (7) 

Here the constraint results from picking the correct pole 
for the residue evaluation in ([5]). It implies that we do not 
get two bound states for every momentum. Surprisingly, 
the branches forbidden by this constraint do not run inside 
the bands or off the unit circle but also in the gap. The re- 
sults are shown in Fig. [3] The overlap of the initial wave 
function with the bound states, i.e. the probability to ob- 
serve a molecule, can be computed from ([5]) and ([7]), see 
Appendix D. It achieves its maximum of | at g = tt. 
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FIG. 3. Dispersion relation cj(p) for the bound state. Parameter: 
interaction phase g in steps of 7r/4. The thick line highlights the 
case g — 7r/4, and dots the case g — tt, for comparison with 
Fig0 

Coin states of the molecule. — If we restrict the Hilbert 
space of the relative motion of the two particles just to 
the bound states, we have a picture resembling very much 
the spectrum of a single walking particle with two internal 
states. Is it actually the same? That is, is there a quantum 
walk with two internal states and only nearest neighbor 
shifts which reproduces exactly the walking dynamics of 
the molecule? For this it is only necessary to compare the 
spectra uo a {p), a = 1, 2, because if these coincide, we can 
find a unitary commuting with translations and mapping 
one walk to the other. Since we have an explicit form of 
the spectrum, we can indeed identify an appropriate coin. 
The walk of the molecules with interaction phase 7 = e %Q 
is unitarily isomorphic to a walk with the second factor in 
([]} replaced by 

1 / 7 V2( 7 -l 
V2( 7 



C = 



(8) 



27 - 1 V V2(7 - 1)7 7 
The exact form of the unitary isomorphism is determined 
by the bound states which are computed in Appendix E. 
The walk uses both branches in ([7]), though only "virtually" 
if one or both of them are forbidden by the constraint. It is 
not clear whether the bound states of general interacting 
quantum walks always allow such an interpretation. 

Velocities. — It is apparent from the right panel of Fig.[T] 
that the molecules are slower than the free particles. The 
square in that figure corresponds to the propagation of one 
site per step. The free Hadamard walk has maximal group 
velocity 1/V2, corresponding to the off-diagonal peaks. 
The group velocity of the molecules (with 7 = —1) is 
duo ± sin p 

-j- = 1 =s=, (9) 

a P v 4 — 4 cos q + sin p 

with the same constraint as in ([7]). The maximal speed (al- 
ways at p = 7r/2) is 1 /3 for the walk with g = n, as shown 
in Fig.[T] For smaller g the maximal speed according to ^ 
approaches one, but becomes forbidden by the constraint. 
Nevertheless it is possible to design molecules, which are 
faster than the free atoms, see Appendix F. 



Prospects for experimental realization. — We plan to re- 
alize walking molecules by extending an existing experi- 
ment [1], in which the walk of single Cs atoms has been 
demonstrated. These happen to be Bosons, so the singlet 
collision case discussed above is not automatically real- 
ized. The internal coin states are taken as the two hyperfine 
states \i) = \F = 3, m F = 3), |t) = \F = 4, m F = 4), 
split by the atomic clock transition at 9.2 GHz. Spin- 
dependent shift operations work, because the up and down 
state see complementary circular polarization components 
of the optical lattice formed by counterpropagating linearly 
polarized waves. Internal coin operations are Rabi rota- 
tions driven by microwaves. 

On-site interactions via s-wave cold collisions can be 
accurately described using pseudopotential methods 03). 
The use of Cs atoms is especially favorable, because its 
large triplet scattering length a T = 2400 a allows fast 
interactions iTTSll . Roughly, given a 3D optical lattice with 
isotropic harmonic confinement uj — 2tt x 30 kHz, the 
state |^o) — (| ti) + I t>l))/v2 acquires a collisional 
phase shift g = tt in about 10 /is. This is just below the 
current step times of our experiment, so could well be in- 
corporated. A further effect of the interaction is a level shift 
for |^o)» which may render the coin operations at 9.2 GHz 
ineffective at collision points, thus suppressing the transi- 
tions to I tt)> I ID- m this wa Y a singlet collision case 
could also be realized for Bosons. For the collision phase to 
be well controlled and coherent, it is necessary to cool the 
motional states to the ground state not only along the axis 
but also in the lateral direction, in which at the moment the 
confinement is much weaker. We plan to enhance the lat- 
eral confinement by an additional blue-detuned doughnut- 
mode laser. Then cooling to average vibrational quantum 
numbers (n) < 6 • 10 -3 in all three directions seems fea- 
sible, which should allow sufficiently many coherent steps 
to clearly demonstrate the molecule binding effect. 
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APPENDIX A: MOLECULES IN HIGHER DIMENSION 

Generally speaking, a point perturbation has less and less 
effect in higher dimension. One therefore expects the for- 
mation of molecules to become more difficult in higher 
dimension. This intuition is basically correct, but since 
points on a lattice have finite thickness, there is no critical 
dimension above which molecule formation disappears. 

From the analysis given in the paper a sufficient con- 
dition is clear: Assume that, for some given total mo- 
mentum p = p , the walk in the difference coordinate 
has a band gap, i.e., there is a z on the unit circle so 
that (W(po,k) — z) is invertible for all k. Choosing 
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r = (I — (see ([6])) we then get a walk with an 

eigenvalue of maximal degeneracy precisely at z. By per- 
turbation theory this will typically split into a collection of 
bound state in a neighborhood of p . This gives a (possibly 
small) range of p and suitable initial preparations, which 
will show the molecule behavior. 

The band gap condition can be met by constructing the 
free walk from a product (not tensor product) of "slow" 
walks in the coordinate directions. Here a slow walk is 
defined as one with nearly constant eigenvalues uj(p). For 
example, consider a walk of the form 

W( Pl , . . . , p s ) = C • ( e ° ipi 6 q 1 )<?!••■ (10) 

■ ■■c ( eiPs \-c 

V e~ iPs / 

By choosing appropriate unitary matrices Ck, a large va- 
riety of s -dimensional walks can be implemented. How- 
ever, when we take all Ck = I, we get a rather trivial 
walk which only has steps along the main diagonal and, 
for all p, only the eigenvalues ±1. Then if we take all 
Ck close to the identity, we get a non-trivial walk in s 
dimensions, whose eigenvalues are . . . ,p s ) ~ 

and cj_(pi, . . . ~ 7T. The same is true for the non- 
interacting pair walk W ® W. Clearly, this now satisfies 
the band gap condition. Intuitively this construction slows 
down the group velocity to such an extent that the pertur- 
bation at a single point can take effect. 

It should be noted, however, that the band gap condi- 
tion is by no means necessary. For example, if we always 
take T time steps together, we are looking at the walk W T , 
whose spectrum, on the one hand, will still show molecule 
states if W does. On the other hand, any bit of continuous 
spectrum will be spread out by a factor T, and eventually 
wraps around the circle, maybe with multiplicity > 1 . The 
molecule states will thus be embedded in the continuum. 

In the continuous-time version of the problem [H the 
band gap condition is always satisfied , since the lattice 
kinetic energy is bounded, and a sufficiently strong point 
interaction can always create a bound state. Hence, as al- 
ready remarked in [8], suitably constructed point interac- 
tions lead to molecule formation even in high dimension. 

APPENDIX B: INTERACTION IN QUANTUM 
CELLULAR AUTOMATA 

Unlike the continuous time case, where one can just add 
Hamiltonian interaction terms, the interaction unitaries in 
the discrete time case have to be multiplied in some or- 
der. This leads to ambiguities for triple collisions. The 
way to treat these is to consider the many-particle system 
as a "lattice gas", a special kind of quantum cellular au- 
tomaton(QCA) lfT6l . This is analogous to a lattice spin 
system with the spin up/down states interpreted as occu- 
pied/unoccupied. The special feature of a lattice gas (as 



opposed to a general QCA) is that there is a notion of parti- 
cle number for each cell, and that the local dynamical rules 
preserve particle number. The characteristic change of per- 
spective is to focus not so much on the different particles, 
but on what can happen in a single cell. 

One standard transition from single particle walk to a 
non-interacting lattice gas is Bose/Fermi second quantiza- 
tion. Since for every basis state \x,a) of the one-particle 
space we have a number operator n Xi(X , we can define the 
number of particles in the cell x as J2 a n x,a- In the Fermi 
case this is clearly bounded by d, but is unbounded in the 
Bose case. This is not quite consistent with the idea of fi- 
nite volume atoms in a potential well, nor with the usual as- 
sumption for QCAs that there should be only finitely many 
states per cell. Indeed the QCA's cell finiteness can be seen 
as a hard core repulsion condition, so a QCA tends to im- 
pose some interaction. There is a procedure for systemat- 
ically introducing an interaction for a walk which is given 
in shift&coin factorization [17]. We illustrate it for the fac- 
torization ([I]) of the Hadamard walk. 

The QCA system consists of a double chain of qubits, 
which we may call the right moving chain and the left mov- 
ing chain. This describes the action of the shift operation: 
In the shift step the entire right (resp. left) moving chain is 
shifted on step to the right (resp. left). The coin operation 
acts at each site. We denote the basis states for all qubits 
by |0) for empty and |1) for occupied, so at each site we 
have the basis vectors \a(3) = |a)i e ft ® |/3) right- m order 
to guarantee particle number conservation, the coin has to 
act in block diagonal form on the number subspaces. The 
empty space 1 00) will just be left invariant. There is no loss 
of generality in this choice, because we are free to choose 
a global phase. On the basis vectors 1 01) and 1 10) of the 
one particle space we just use the coin operation of the free 
walk. Then it is easy to see that the subspace of total par- 
ticle number 1 is invariant, and isomorphic to the free one- 
particle walk. The only free parameter in this construction 
is the phase 7 for the two-particle state |11). Again, it is 
easy to check that on the 2-particle subspace we exactly get 
the singlet collision case discussed in the paper. However, 
the construction also extends to arbitrarily many particles. 

This construction can be extended to other walks with an 
explicit decomposition into shift and coin operations (al- 
lowing also for many factors). Since such decompositions 
are not unique there are some trade-offs possible[17j: For 
example, one can use a decomposition in which only one if 
the d states is ever shifted. The multiple chain picture then 
only has a qubit chain for transport, and a stationary chain, 
with an empty state and, for particle number 1, the remain- 
ing (d — 1) states of the one-particle walk. This leads to 
a QCA with 2d states per cell. However, the interaction 
neighborhood size of this QCA may be much larger than 
that of the one-particle walk. In contrast, taking a moving 
chain for every one of the d states suggests a picture with 
2 d states per cell, and in this case (at least in ID) the neigh- 
borhood can be arranged to be the same as for the walk. 
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Further options are introduced by allowing more than sin- 
gle occupation of the one-particle coin states, as in the free 
Bose walk. 

Which of these many possibilities is realized is, of 
course, a question that has to be answered by analyzing 
the physics of the given system. The point to be made here 
is only that the language of QCAs makes clear what physi- 
cal information is required to specify an "interaction" com- 
pletely. 



Then 



U = I - R(z) 



-I 
-I 



i 



\{l + ik) 
-l + ik 



E K (dk) 



l + ik 



E K {dk). 



Now, since | — 1 + ik \ = | 1 + ik | the integrand has absolute 
value 1, and so the integral represents a unitary operator. 



APPENDIX C: A UNITARITY LEMMA 



APPENDIX D: OVERLAP WITH BOUND STATES 



In this appendix we prove the claim that the operator 
(1 — R(z)~ l ) appearing in ([6]) is unitary. We will do this 
by proving the following, more general claim: 

Consider a unitary operator W on a Hilbert space % and 
a point z on the unit circle (z G C, \z\ = 1) not in the spec- 
trum of W, which just means that the inverse (W — 
exists and is a bounded operator. Now let N be an arbitrary 
projection and set R(z) = N(W — zl)~ 1 WN, consid- 
ered as an operator on NT-L. We claim that this operator is 
invertible and that U = 1 — i?(z) _1 is unitary. 

The crucial step for the proof is to show the following 
relation 



R(z) + R(zY 



(11) 



We can show this for the special case N = I, and 
then multiply the above equation with N from both sides. 
Since W is unitary it has a spectral decomposition W = 
JQ E w {dVt), where E w is the spectral measure on the 
unit circle. We can then express the operator ( [TT] ) directly 
in the functional calculus (we use f for hermitian adjoint 
and an overbar for complex conjugate): 



R(z) + 



+ = I E w (dft) 



z) - z 



I \n-z\i 

2 - 2U{Tlz 



Ew(<m) 



I 



2 - 2^(ftz) 



E w (dQ) 



This determines the hermitian part of R(z), and we can 
write R{z) — |(1 + iK) for some hermitian operator K 
on NU. 

For the next we use the functional calculus of K = 
J k E K (dk), with the spectral measure of K, which is now 
supported on the real axis. We note that due to the projec- 
tion N, which need not commute with W, this measure 
cannot be easily obtained from the spectral measure E w . 



The time evolution of a two-particle quantum state under 
the considered quantum walk W is essentially determined 
from its overlap with the bound states ^(p) in Fourier 
space. This is also important for experimental implemen- 
tations since it describes how efficient some preparation, 
e.g., "both particles in singlet state at the origin", is for 
generating molecules. In other words, we want to estimate 
the splitting in Fig. [T] between the total probability near the 
diagonal and the rest, i.e., the probability for ballistically 
moving away from the diagonal. Since the total momen- 
tum p is a conserved quantity we omit the p-dependence of 
W in the following discussion. 

First, let us calculate the normalization of the eigenstate 
^ of the quantum walk corresponding to eigenvalue z as 
determined from ([4]). Abbreviating tpr — (1 — T)^we get 



(Vt| -z 



d s k 



"Vr|| 2 
W(k) 



(2tt) s (W(k) - zf 



= (*!* 



d 



(12) 



where we used the definition of R(z) from equation ([5]) 
and the unitarity of W(k). This yields now the probability 
to capture the walker initially prepared in an internal state 
$ located at the origin in the bound state 



P 

± cap 



^\R(z)^ 



w \\w ' (13) 

since ij) respects ([6]). We want to evaluate fl3] ) for our stan- 
dard example, the Hadamard Walk with singlet collisions 
and the singlet state ^_ as the initial state, where ^ r is 
given by (1 — 7) By ( [T2| ) we have to calculate the re- 
ciprocal of 

m * = _|i _ ^ Z j-^_\R{ Z )^_) = -iT 1 ^ 

where we used ([6]) and the convention z — e lUJ in the last 
step. Defining 77 via e l7] = — - one can deduce from ([7]) 
the relation 



cos 77 = cosp — 2 cos oj 



(14) 
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Hence, we get for the capture probability P cap j 



in the bound 

state corresponding to the eigenvalue z the expression 

-l 



1 + 2- 



smo; 



smr] 



(15) 



Note that this is a function of the total momentum p since 
r] as well as uj depend on it. Analyzing this dependency, 
we see that as uj approaches a band edge, sin r] tends to 
zero and therefore such p values are suppressed as Fourier 
components of the bound state. That behavior can also be 
seen in Fig. |4j 




FIG. 4. The figure shows the capture probability P ca p for the cj+ 
branch (left panel) and the sum of both branches (right panel) as 
a function of the total momentum p and of the interaction phase 



For fixed interaction phase 7 = e %9 the ratio between the 
occurrence of molecules and unbound two particle states in 
an experiment is determined from the integral of P cap with 
respect to the total momentum p. As can be seen from Fig. 
[5] there is a fairly high probability to observe molecules 
with a maximum of 2/3 at g — tt. 




FIG. 5. The integrated capture probability P ca p is plotted over 
the interaction phase 7 
at g = 7r. 



e %g . It reaches its maximal value of 2 /3 



Expression ( [T5] ) also allows us to study the long time 
behavior of the molecule's position distribution scaled by 
1/t, where t is the discrete time parameter. As outlined 
in corollary 7 in |[T2l . this asymptotic distribution can be 
evaluated with knowledge of the group velocity ([9]) and the 
capture probability ( [T5] ). The result is shown in Fig. [6] 

APPENDIX E: ANALYTIC FORM OF BOUND STATES 

The goal of this section is to derive an explicit formula 
for the bound states of a quantum walk with a point defect 




FIG. 6. The left panel shows the molecule's asymptotic position 
distribution depending on the interaction phase 7 = e zg . The 
yellow line corresponds to an interaction phase of g = 0.8, the 
red line to an interaction phase of g = tt. The flat distribution 
profile for g — 0.8 matches the flat position distribution of the 
bound state shown in the right panel for t — 50 time steps of the 
interacting quantum walk. The two peaked caustic behavior in 
the asymptotic position distribution for g — tt can be compared 
to the corresponding panel in Fig. [T] 



at the origin. We then apply this procedure to our stan- 
dard example, the Hadamard walk with singlet collisions. 
This allows us to derive an explicit form of the p dependent 
transformation needed to construct the virtual walk ^ of 
the molecule as explained in the main part of the paper. 

The components ^ x of the unnormalized eigenvector \I/ 
corresponding to lattice site x can be determined from ([4]), 
with the substitution v — e %k the formula reads 

**=hS^ {w{v) - z) ~ lw{v) ^ 

= 2^/^( I + ^)-^)" 1 )^- (16) 

Here, we have to substitute the correct form ([7]) of the 
eigenvalue z — e luJ . The initial state in our example is 
the singlet state which effectively reduces the problem 
to the Fermi sector. Of course, the Fermi symmetry is re- 
flected in the eigenstate If we denote the operator which 
exchanges tensor factors by F, that is, F77 ® x — X ® V 
for all vectors 77, x E C 2 , then Flf(v)F = W^' 1 ) and 
hence ^ x = -F^_^ for all x e Z. 

In order to determine ^ x for negative x we use resid- 
ual calculus which requires knowledge of the poles of the 
integral kernel in ( [T6] ). The operator (W(v) — z)~ x may 
have singularities at the zeros of the quadratic polynomial 
p(v) = v det(W(v) — z) and possibly at v = 0. Conjuga- 
tion of W(v) —z by F yields the equality det (W(v) — z) — 
&et(W(v~ l ) — z), which implies that the singularities of 
(W — z)~ x are either of first order and inverses with respect 
to each other or p(v) is constant, in which case the operator 
(W — z)~ x has a pole of first order at v = 0. Analyzing 
the case where p(v) is not constant we get the following 
expression for the singularity with modulus smaller than 
one, denoted by Vi, 

Vi — — cos 77 — cot - sin 77 
provided sin(?7 + g) • sin(?7) < , (17) 
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with 77 as in ( [14] ). Clearly, ^ — V>- an d if x < only 
2(W(?;) — z) contributes to the integral ([16]). Since the 
integral kernel is analytic at v = it remains to calculate 
the component wise residual operator 

R V1 = Res (v^ziWiv) - z)' 1 ) 

v—tvi 

at V\. The coefficients of the normalized eigenstate \& in 
terms of this operator read 

( (l-j)v[ xl R Vl xl;_ , x<0 
&x = Pcap • < i>- , a = . 

I -(l- 7 )^Fi^_ , x>0 

Note that since the modulus of V\ is smaller than one the 
state \I/ decays exponentially in \x\. This implies that the 
distance of the two particles in the molecule state is expo- 
nentially concentrated around zero. 

In the case that the singularities of (W(v) — z)~ x are 
at Vi = the eigenstate \I/ is strictly localized on a finite 
set of lattice sites. This is because for Vi = the integral 
kernel in fT6] ) is analytic if \x\ > 1, hence *f? x = for 
these cases. For our example this happens exactly at p = 
±uj(p). This property leads to an interesting feature of the 
considered quantum walk: By engineering the initial state 
of the quantum walk on few lattice sites it is possible to 
generate states which are close to a momentum eigenstate 
at points p = ±ou(p). Here, the spread of the initial state 
in the relative coordinate is bounded by the spread of the 
eigenstate 4? x and in the center of mass direction by the 
desired accuracy of the momentum preparation. 

With the help of the eigenstates ^± corresponding to 
the two branches uj± we can construct the one dimen- 
sional quantum walk mimicking the time evolution of the 
molecule explicitly. This is done by identifying the states 
with the p dependent eigenvectors of some one dimen- 
sional quantum walk with the same dispersion relation oj±. 



APPENDIX F: FAST MOLECULES 

So far only the two particle Hadamard Walk with singlet 
collisions was considered as explicit exemplification of the 
ideas of this paper. From Fig. [I] and Fig. [6] one could think 
that the bound states are slower or at most as fast as the free 
walkers. 

In this section we will give two numerical examples to 
show that this is not true in general and that it is possible 
to design interactions that allow for molecules spreading 
faster then the free walk. Since in the case of singlet colli- 
sions the only freedom lies in the selection of the interac- 
tion phase 7 we turn to the Bose case where we can choose 
an arbitrary three by three unitary matrix in symmetric sub- 
space. 

A trivial possibility to generate such fast molecules is 
choose an interaction matrix that counteracts the coin op- 
eration. In that case the walk operator for the two particle 



Hadamard walk is given by 

W r = (Wi ® Wi) ((I -N) + H® HN) 
= (S ® S)(H <g> H{1 — N) + IN) , 
which acts on the walkers only with the shift at the collision 
point. If we now prepare the particles at the same lattice 
side for example in the state | ft) they W1 U be just shifted 
to the left with maximal velocity. 




4-50 



FIG. 7. The two panels show two particle Hadamard walks with 
different bosonic interactions. In both cases the molecule state 
extends beyond the green square indicating the maximal velocity 
of the free Hadamard walk. 

Two nontrivial examples of fast molecules are given in 
Fig. [7] were we have choosen two unitary interactions for 
which the bound state in the symmetric subspace spreads 
faster than the individual particles of the free quantum 
walk. Moreover it is even possible to design the interaction 
in such a way that the resulting molecule travels nearly at 
the maximal possible velocity as in the trivial case. 



[1] M. Karski, L. Forster, J.-M. Choi, A. Steffen, W. Alt, 

D. Meschede, and A. Widera, Science 325, 174 (2009) 
[2] H. Schmitz, R. Matjeschk, C. Schneider, J. Glueckert, 



M. Enderlein, T. Huber, and T. Schaetz, Phys. Rev. Lett. 
103, 09 0504 (2009) 
[3] F. Zahringer, G. Kirchmair, R. Gerritsma, E. Solano, 
R. Blat t, and C. F. Roos, Phys. Rev. Lett. 104, 100303] 
(2010)| 



[4] M. A. Broome, A. Fedrizzi, B. P. Lanyon, I. Kassal, 

A. Aspuru-Guzik, and A. G. White, |Phys. Rev. Lett. 104,| 
"3602 (2010)] 
[5] A. Peruzzo, M. Lobino, J. C. F. Matthews, N. Matsuda, 

A. Politi, K. Poulios, X.-Q. Zhou, Y. Lahini, N. Ismail, 

K. Worhoff, Y. Bromberg, Y. Silberberg, M. G. Thompson, 

and J. L. OBrien, Science 329, 1500 (201 0)| 
[6] Y. Bromberg, Y Lahini, R. Morandotti, and Y Silberberg, 

Phys. Rev. Lett. 102, 253904 (2009) 
[7] A. Schreiber, K. N. Cassemiro, V. Potocek, A. Gabris, P. J. 

Mosley, E. Andersson, I. Jex, and C. Silberhorn, Phys. Rev. 

Lett. 104, 050502 (2010)1 
[8] K. Winkler, G. Thalhammer, F. Lang, R. Grimm, J. H. Den- 

schlag, A. J. Daley, A. Kantian, H. P. Biichler, and P. Zoller, 

Nature 44 1, 853 (2006)| 



9 



[9] M. Stefanak, S. M. Barnett, B. Kollar, T. Kiss, and I. Jex, 
|New J. Phys. 13, 033029 (20TT)1 
[10] N. Inu i, Y. Konishi, andN. Konno, [Phys; Rev. A 69, 052323) 
(2004)| 

[11] Y. Omar, N. Paunkovic, L. Sheridan, and S. Bose, Phys. 

Rev. A 7 4, 04230 412006)1 
[12] A. Ahlbrecht, H. Vogts, A . Werner, and R. Werner, J. Math. 

Phys. 52, 042201 (20TT)1 



[13] T. Kato, Perturbation theory for linear operators, 2nd ed. 
(Springer, 1976). 



[14] R. Stock and I. H. Deutsch, Phys. Rev. A 73, 032701 (2006) 



[15] P . J. Leo, C. J. W illiams, and P. S. Julienne, Phys. Rev. Lett.| 
85, 2721 (2000)| 

[16] B. Schumacher and R. F. Werner, "Reversible quantum cel- 
lular automata," arXiv:quant-ph/0405174. 

[17] H. Vogts, "Discrete time quantum lattice systems," Doctoral 
Thesis, Technische Universitat Braunschweig (2009). 



